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Abstract 


We are concerned with constrained Chebyshev approximation problems of the type 


min 

pen n : P (c)=i 


max IK*) I 

zes 


( p ) 


Here H n denotes the set of all complex polynomials of degree at most n, £ is any ellipse 
in the complex plane, and c £ (F \ £. Such approximation problems arise in the context 
of optimizing semi-iterative methods for the solution of large, sparse systems of linear 
equations Ax — b with complex non-Hermitian coefficient matrices A . The problem of 
obtaining optimal ploynomial preconditioners for conjugate gradient type methods for 
Ax = b also leads to problems of type (P). In this paper, we introduce a new family of 
polynomials q n {z\ c), n 6 iV, c £ <F\£, as the polynomials which are optimal for a modified 
version of (P) with n n replaced by a certain subfamily. Some simple properties of q n are 
also listed. We then derive a necessary and sufficient condition for q n to be the extremal 
polynomial of the approximation problem (P). Finally, it is shown that q n is indeed optimal 
for (P) for all fixed n whenever the distance between c and £ is sufficiently large. Results 
of some numerical tests are presented. 


1. Introduction 

In this paper, we are concerned with constrained Chebyshev approximation problems 
of the type 

(D n (r, c) :=) min max |p(*)| . (1) 

pen n ;p(c)=i ze£ r 

Here H n denotes the set of all complex polynomials of degree at most n, 

S T := {z e <T I \z - 1| + \z + 1| < r + -} , r>l , (2) 

r 

is any ellipse (including its interior) in the complex plane with foci at ±1, and it is always 
assumed that is any ellipse (including its interior) in the complex plane c 6 <F\ £ r . Since 
all polynomials p 6 H n with p(c) = 1 can be parametrized in the form 

p(z) = l-(~n(z-c) + 7a(z - c) 2 + f7n (*- c )") > 71,72,. •• ,7n € <T , (3) 

the condition c G <F\ £ r guarantees that Haar’s condition is satisfied. Thus, there always 
exists a unique optimal polynomial for (1) which will be denoted by p n (z\ r,c) in the 
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sequel. However, these extremal polynomials are explicitly known only for special cases. 
The solution of (1) is classical for real c: 


Pn{z\ r,c) 


Tnjz) 

r»(c) 


cG M\£ r , 


(4) 


where T n is the nth Chebyshev polynomial (of the first kind). 

Constrained approximation problems (1) with complex c arise in the context of opti- 
mizing semi-iterative methods for the solution of non-Hermitian systems of linear equations 
(e.g. Manteuffel [4] and Eiermann, Niethammer, and Varga [1]). Mainly motivated by this 
application, in some recent papers, problem (1) was studied for complex c and the opti- 
mal polynomials were found for certain special cases. For n = 1, Opfer and Schober [6] 
obtained a complete solution of a more general version of (1) with £ C <F any compact set 
not containing c. For ellipses, their result can be rewritten in the following form: 


pi(z; r,c) 


Bz + i sin 7 
A (B cos 7 + i A sin 7) 


(5) 


where 


c 

with 0 < 7 < 2n and 


A cos 7 + 1 B sin 7 (G d£ R ) 


( 6 ) 


' 4 - K* + s) ’ b = ’ R>r zi , 

(by d£ R we denote the boundary of S R ). Freund and Ruscheweyh [3] investigated (1) for 
the case r = 1 of the line segment £ x = [-1, 1], They determined p 2 (z; l,c) for arbitrary c 
and Pn(z; 1, c) for n G N and purely imaginary c. In both cases, the optimal polynomials 
are suitable linear combinations of T n , T„_ x , and T n _ 2 . Finally, Fischer [2] showed that 
for nondegenerate ellipses £ r , r > 1, and purely imaginary c the normalized Chebyshev 
polynomial (4) is optimal for (1), if n is even and |c| is sufficiently large compared to r. 

Note that, except for the cases solved in [3], all the other explicitly known optimal 
polynomials are of the form 


9(2) 


T n (z) + a 

T„(c) + a 


a € <F . 


(7) 


It is thus natural to ask, whether polynomials of type (7) lead to explicit solutions of (1) 
also for the case of general complex c and n G IN . The purpose of this note is to answer 
this question. 
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The paper is organized as follows. In Sect. 2, we introduce a new family of polynomials 
q n (z; c), n E iV, c € (P \£ r , as the polynomials of the form (7) with minimal uniform norm 
on £ r . Some simple properties of q n are also listed. In Sect. 3, we derive a necessary and 
sufficient condition for q n to be the extremal polynomial of the approximation problem (1). 
Finally, Sect. 4 contains the main result of this paper. We show that indeed p n (z: r, c) = 
q n (z ; c) for all fixed n E W, r > 1 and all c E <F whose parameter R in the representation 
(6) is sufficiently large, i.e. R > i2o(n,r). An explicit formula for i?o( n ? r ) is given. 

2. A Class of Extremal Polynomials 


Throughout this paper, let n E JV, r > 1 ° r be the ellipse defined in (2), and it is 
assumed that eg W\£ r with representation (6). We will make use of the parametrization 

z r (<f>) = a cos<f> + i b sin<£ , <f> E IR , 

of the boundary d£ r of £ r . Here a := ax, b := where 

<■*»;(<•* + -X) and ‘* : = 5( r ‘ " 4) ’ * = 1-2. • (8) 

2 r K 2 r K 

Tk(z) denotes the kth Chebyshev polynomial which by means of the Joukowsky map is 
given by 

Tk(z) = V + -L) , z = i(r + i) . (9) 

2 v K 2 v 

By (6) and (9), one has 

Ck := Tfc(c) = Akcos(ky) + i B k sin(kj) , A: = 1 , 2 ,... , (10) 

where 

A * ■= \( Rk + -k') and ■ (in 

The relations 

A? k — B\ = 1 , a\ — b\ = 1 , k = 1, 2, ... , 

will be used repeatedly in the sequel. Moreover, note that, since R > r, 


Ak > ak , Bk> bk , 

We consider the extremal problem 

(M n (r,c) :=) min max 


& — 1,2,... 
r w (z) + a 

Tn(c) + Or 


( 12 ) 
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Since w = T n (z ) maps £ r onto £ r n , (12) is equivalent to 


min max 

p€lli :jj(c n ) = l w££ r n 


|p(«;)| . 


Thus, by (5) and (10) (for k = n ), 


Qn(z]c) := pi(T n (z); r",c n ) 


B n T n {z ) + i sin(n7) 
A n (J3 n cos(n7) -f z sin(n7)) 


( 13 ) 


is the unique extremal polynomial of (12). Next, we determine M n (r,c ) and the corre- 
sponding extremal points, i.e. :gf r with 


kn(z’,c)\ = M„(r,c) . 

From the maximum modulus principle it follows that all such points he on d£ r . By (9) 
and (8) (both for k = n), one has 


T n (z r {4>)) = cos(n<£) + i b n sin(n^) . 

Using this identity, we deduce from (13) the relation 

i i t ±\ m2 a n ( , {B n sin (n<f>) — b n sin(« 7 )) 2 \ 

\<lTx{Zr{<t>)-,c )\ 2 = -77 ( 1 - 2 /m l - 2/ U ’ ( 14 ) 

A n \ a n(B^ + snr(n 7 )) ) ’ 

Therefore M n (r, c) = a n / A n , and the extremal points are just the z r (<j>) with <j) satisfying 


We set 


and define V’n hy 


Note that 


B n sin (n<t>) = b n sin(nj ) . 

d n := -^-sin(n7) 

o n 


• I J 7T 7T 

SinV’n = d n , -- < V’n < - - 


M.I < w < 1 

All solutions (mod 2t>) of (15) are then given by 


(15) 

(16) 

(17) 

(18) 


<tn = -* + (-1)'— , l = 1,2, . . . ,2n . 

n n 

Remark that for r > 1 (resp. r = 1) this leads to precisely 2n (resp. n + 1) distinct 
extremal points of q n on d £ r . We summarize these results in the following 
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Theorem 1. q n (z\ c) is the unique extremal polynomial of (12), and the corresponding 
minimal norm is 


M n (r, c) = 


r n + 1 / r” 


R n + 1 / R n 

On S r , r > 1, q n (z;c ) has precisely 2 n extremal points: 

z i = ~(r + -) cos <f>i + -(r --) sin <f>i , <fn = -n + (-1)' — 
It It n n 

The extremal points of q n (z ; c) on Si = [—1, 1] are 

In 

z\ = cos — , 1 = 0, . 

n 


l = 1,2,. ..,2 n 


Remark l.The optimal polynomial of (12) is identical for all S r , 1 < r < R. M n (r,c ) 
depends only on the parameter R of Sr, but not on the position of c on 8Sr. 

The family of polynomials q n (z; c ) also leads to upper and lower bounds for the min- 
imal deviation D n (r,c ) of (1). 

Theorem 2. Let r > 1, c E 8 Sr , R > r. Then, 


Dn(r,c) < 


r n + 1 / r” 

fl n + 1 / il” 


n — 1,2,... 


and 


for all n satisfying 


D n (r,c) > 




(Bn + bn |sin(n 7 )|) 2 
a U B l + sin 2 (n7)) 


|sin(i» 7 )| < b n B n . 


(19) 


( 20 ) 

( 21 ) 


Remark 2. Clearly, (21) is true if n is sufficiently large. 

Proof. (19) is an immediate consequence of Theorem 1. A standard technique (Trefethen 
[8], ManteufFel [4]) to obtain lower bounds for complex approximation problems is based 
on Rouche’s theorem. Applied to (1) and q n , this yields 


D n (r,c) > min \q n (z;c)\ 

z£dt r 


( 22 ) 


if it is guaranteed that all zeros of q n are contained in S r . In view of (14), the right-hand 
side of (22) is just the bound stated in (20). By (13), the zeros of q n are the solutions of 
the equation 


Tn(z) = 


— X 


. sin(n7) 
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Using (9) (for k = n), one easily verifies that all these solutions he on the boundary d£ p 
of an ellipse of type (2) whose parameter p > 1 is defined by 



|sin(n 7 )| 


B 


n 


Therefore, d£ p (and hence the zeros of q n ) is contained in £ r iff /?„ < b n . This concludes 
the proof of Theorem 1. ■ 


3. A Criterion for Optimality 

As mentioned in the introduction, it is known that 

p n (z-,r,c) = q n (z;c) , z £ <T , (23) 

for some special cases as n = 1 or c £ 1R \ £ r . In this section, we present a necessary and 

sufficient condition for (23) for the general case n £ N, c £ <F\ £ r . This criterion allows 
to check (23) by computing 2 n real numbers for which explicit formulas are derived. 

First, consider the case r = 1 of the degenerate ellipse £\ — [ — 1, 1], It was shown in 
[3] that p n (z; 1, c) has precisely n + 1 extremal points 

1 = z 0 > Zi > • • • > z n = —1 
and there is a s n £ (T such that 

p„(z/;l,c) = s„(— 1)' t ~ — “7 , / = 0,1, . . . ,n . (24) 

\zi — c I 

By Theorem 1 and (13) (with z = z/), q n (z'c) has the extremal points 

*i = cos ~ and <ln(zi;c ) = t n ((-l) 1 + i — n ( n7 ) ) , / = 0, 1, . . . , n , (25) 

for some i n £ <T. By comparing (24) and (25), it is straightforward to verify that, for r = 1, 

(23) holds iff n = 1 or c G 1R \ £ r - So, except for the already known cases, q n {z\c) is not 
optimal for (1) with r = 1. 

Therefore, for the rest of this paper, we assume that r > 1. By Theorem 1, the 
extremal points of q n (z', c) on £ r axe 

l ib 

Zi := a cos <f>i + i b sin 4>i , <f>, := -x + (-1)— , / = l,2,...,2n , (26) 

n n 

with V’n defined by (17) and (16). We list some properties of the points (26), which will 
be needed for the derivation of the main result of this section, in the following 
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Lemma 1. a ) For / = 1,2,..., 2 n: 


sin (nfa) = <f n , cos (nfa) = (-l)'^/! - , 


and 


where 


q n (z,;c) = < w ((-l)* +» -- q "_l n ) 

6 nV /l - d* n 

a n B n yJ 1 - <£2 


(27) 


(28) 


<n = 


bj For j = 0,1,..., 2n: 


2n 


i=i 


A n (B n cos(n 7 ) + i A n sin(n 7 )) 


r i it j = o 

i d n if j = n 

l-2d 2 n if j = 2n 

0 otherwise 


=2n x 


(29a) 


and 


2 n 




f=l 


if j = n 


2 id n if j = 2n 
\ 0 otherwise 


( 296 ) 


Proof. (27) follows immediately from (17) and the definition of 4>[ in (26). (28) is obtained 
from (13) (with z = z\) by using (9) (for k = r^), (16), and (27). 

We now turn to the proof of part b). Recall that 


. t. f n if j e n2Z 
to if j £ n2Z 


(30) 


Let 0 < j < 2n and 8 = ±1. Since, by (26), 


t-hi + 'k i( , = 2k 


<f>i = 


71 71 


(2k - 1 ) 7T V’n 


if / = 2k - 1 


71 


71 


and with (30), we get 


1=1 


k = 1 


k=l 
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, ij±IL c - **; _ , , JtlL Jt 

== (e J n +<5e » e ' ?n )y e *» = n 


ii —if x 


X < 


*= 1 


f 1 + 8 if j = 0 

e*^n _ &,-**„ if j = n 

e 2itl> n + ^ e -2.>„ if j _ 2n 

0 otherwise 


Using (17), one easily verifies that these are just the formulas (29a) (6 = 1) and (29b) 
(6 = - 1 ). ■ 

In view of (3), (1) is a linear Chebyshev approximation problem: We seek the best 
uniform approximation to f(z) = 1 on £ r out of all functions of the linear space 

n„(c) := {p e n„ | p(c) = o} . 

Therefore, the characterization of best approximations due to Rivlin and Shapiro [7] can 
be applied. The following criterion results: 

Criterion 1: q n (z-,c) is the optimal polynomial for (1) iff there exist nonnegative real 
numbers a\,<T 2 , ... ^cr^n (not all zero) such that 

2n 

^2 <t, q n (zp, c) p(z,) = 0 for all p G D„(c) . (31) 

f=i 

We now determine all real which fulfill (31). Note that q n (zp,c ) is given 

explicitly in (28). Furthermore, II n (c) is spanned by the polyno mials 


and, by (26), (9), 


Tk(z) — ck , fc = l,2, ...,n 


T t (r,) = i(r fc e i ^' + ^e-*^') 

Z r K 


Thus, (31) can be rewritten in the form 


2 n 


£>,((-!)' -«*)(rV**+r *e - 2c») - ° , k = 1, 2, . . . ,n , 


( 31 ') 


/=i 


where 


e n * — 


a n d n 


Next, we remark that any numbers <7i, . . . , <r 2n G JR admit a representation of the type 


(32) 


cos (j<h) + sin (j<j>,)) 


i=o 

n 


(33) 


= !>>«"'*+* • ' = 1.2.-.2n , 

i = o 
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with real numbers A j, fij, j = 0, . . . , n, fio •= 0, and 



This follows from the fact that the linear space spanned by 


(34) 


1 , cos<^> , cos(20) , ... , cos (n<f>) , sin</> , sin(2<£) , ... , sin(n^) 

satisfies Haar’s condition on any interval of the form [a, a + 2i r), a £ iR, and since, by (26) 
and (17), the numbers <f>t y l = 1, . . . , 2n, are distinct and all contained in such an interval. 
By (33), (31') leads to a system of equations for i^o, . . . , : 

n 2 n 

+P7c _i,>, )(r*c i *** +r-*e-'**‘ -2c*) = 0 , 

j — 0 /=1 

fc = 1, 2, . . . , n . (31") 

A routine calculation, making use of (29a, b), (32), and (34), shows that (31") reduces to 

u n - k r k (b n + + P^ri-^-(6 n - r n d 2 n ) - ia n d n ( ^ + F*r*) 

— 2c*(A n 6 n (l d n ) iund 7 i(Ao ~b d n /x n )) , fc = 1, 2, . . . , n 1 , (35u) 

and, for k = n, to 

^n(^n “b ic n d n )(Ao ~b d n // n ) (1 d n )(6 n c n -b id n )A n = 0 . (356) 

Note that Aq and /i„ only occur in the combination 


T — a n (Ao ~b dnl^n) 


(36) 


moreover, we set 

A := (1 — d 2 n )\ n . (37) 


By taking its real and imaginary part, respectively, each of the complex equations (35) 
yields two real equations. Using (34), (8), (36), and (37), we thus arrive at 

(<zjt6 n b n — kd n )A n _* -|- d n bk(i n fj>k = 2(.Re c & )6fj, A -f" 2 (7m )d n t , 

(35'a) 


d] c (ind n \} i “I - (6^6 n “b o>n— k k — 2(7?e Cfc ^d^T 2(7mc^)6 n A , 
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for k = 1, . . . ,n — 1 , and 

(b n - (Imc n )d n )T - (Re c n )b n \ - 0 , 

(Re c n )d n T - (b n (Im c n ) + d n ) A = 0 . 

With (16) and (10) (for k = n), the two equations of (35'b) can be written as 

cos(n 7 )(r cos(n 7 ) - A A„) = 0 , sin(n 7 )(-r cos(n 7 ) - AA n ) = 0 . 

Therefore, the 2x2 system (35'b) is of rank 1 and its solutions are described by 

t e 1R . 


(35'6) 


cos(n7) 

A = ; T 




(38) 


Now assume that r € JR is arbitrary, but fixed, and let A be defined by (38). It remains 
to solve the system (35'a) of 2(n - 1) linear equations for the 2(n - 1) unknowns A* and 
fik, k = 1, . . . , n - 1. First, we note that, by combining the first equation of (35'a) with 
the second one of (35'a) (with k replaced by n - k), the system (35'a) is equivalent to the 
n — 1 decoupled 2x2 systems 


where 


and 


C k 


= 26 “ T u) 

, k = l 

a kbn — 

dn bfcdj 


(39) 


b n -kb n + a k (P n 


0 

n / 


h = (ReCk) ^pl +(ImCk) ^jl , )k=(Rec ^ k fJ^ll_ (ImCn _ t) M 


A n 


B n 


An 


(4 °) 

Here, the formulas (40) were obtained by using (38) and (16). With (8), it is easily verified 
that 

det C k = a k b n -k(b^ + d^)(l — d\) . (41) 

Thus, in view of (18), all matrices C k , k = 1 , . . . , n — 1 , are nonsingular, and by Cramer’s 
rule we deduce from (39) and (41) that 


where 


A* = rA*(l) , fi k = Tfi k (l) , (42) 

( ( b n , d?n\t 0,nd n b n - k \ , 

(> & + - <) 4 * a*-* 9 ”-V (43a) 
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and 


M 1 ) = 


2b n 


/ b n d n a n d n fin—* , 

( t _ )Qk “ 7 Jk 


(436) 


(62 + dl)(l - <P n ) V v 6„_ fc a k a k 6 n _fc 

& = 1, . . . , n — 1. Finally, note that, by (27), (36), (37), and (38), summing up of the first 
(j = 0) and the last (j = n) term in (33) yields 

A 0 + An cos(n<f>i) + Hn sin (n<j>i) = r(— + _c os(n 7 ) ^ = (44) 

a n y 1 — d* 

Summarizing, we have proved that the set of all solutions a — (a 1 , <r 2 , . . . , <r 2n ) T € iR 2n of 
(31) is given by the one dimensional linear space 


a = — a* , r € M , 


Or 


where, by (33), (42), and (44), 


a t : = 1 + ("I)* 7==f 'J L + a " X](A*(l)cos(fc<£{) +/ijfc(l)sin(^/)) , 

V 1 - d n Afl fc=1 


/ = l,2,...,2n . 


(45) 


Hence, Criterion 1 can be restated as follows. 

Theorem 3. Let n £ iV, r > 1, c £ <F\ £ r . Then, the polynomial (13) q n { z \ c) is optimal 
for (1) iff the numbers (45) <j*, / = 1,2,..., 2n, are either all nonnegative or all nonpositive. 

Remarks 3. For given n, r, c £ 8£r, R > r, the numbers erf, 7 = 1,... ,2n, can easily 
be computed numerically by means of the formulas (6), (8), (11), (16), (26), and (43). We 
have done that in a number of cases. These numerical tests indicated that the polynomials 
q n (z]c) are indeed optimal for (1) whenever R (for fixed r, n) resp. n (for fixed r,R) is 
sufficiently large. We were not able to characterize explicitly all n,r,R for which q n is 
optimal. However, in the next section, a necessary condition for the optimality of q n is 
derived. 

4. For the simplest case n = 1, the sum in (45) does not occur. It is easily verified that 
R > r guarantees *r > o. 7 = 1,2, and thus we have reobtained the result of Opfer and 
Schober [6] for the case n = 1. 

5. It follows from Meinardus’s invariance theorem [5, Theorem 27] that the extremal 
polynomials of (1) corresponding to c and its reflections c resp. — c on the real resp. 
imaginary axis are connected through 


p n (z;r,c) = p n (z;r,c) resp. p„(z;r, -c ) = p n (-z;r,c) , z G (F . 
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This symmetry is also reflected in the following relations for the numbers (45). For fixed 
n and r, we consider of = of(c) as a function of c. Then, 

*t( s ) = <rL-i( c ) > / = 0,l,...,2n , 

and 

\ <-l( c ) > / = 0,l,...,n 

1 a 3 n-l( c ) > / = n + l,...,2n 

where of := of n . These identities can be verified by a routine calculation using the 
definition of of. 


4. Optimal Polynomials for the Constrained Chebyshev Problem 


In this section, we present a simple inequality involving n, r, R which guarantees the 
optimality of q n for (1). For that purpose, a lower bound for the numbers (45) is derived 
which finally leads to a necessary condition for the nonnegativity of of , / = 1,2,... ,2n. 

Throughout this section, it is assumed that n > 2 , R > r > 1 , and that c € 8 £r is 
represented in the form (6). Moreover, we recall the definitions of a*, 5 * (in (8)), Ak, Bk 
(in ( 11 )), d n (in ( 16 )), and /*, <7* (in ( 40 )). In the following lemma, some estimates, which 
will be used in the sequel, for these numbers are listed. 


Lemma 2. a) For k = 1, 2, . . . , n — 1: 


b) 


i/jt I < M 


£n(l ~ d 2 J 

A 2 n~ a2 n 


15*1 < A, 


* 


g.(i - 4 ) 

Al-al 



< 


4 r 5 


(?r 


< 


(r 4 — l)(f? — r) v r 
4 r 4 

(r 4 — l)(Hr — 1) 


R n 


B 


k=l 


k 

^k^n—k 


+ 


bn-k 

Qn—kbk 


)Ak 


< 


2r(2r 2 + 1) R 
(r 2 -l)(J2-r)^r ’ 


( 46 ) 


Proof, a) By Cauchy’s inequality, it follows from ( 40 ) that 


\fk\ < \ck\y/g{x) , I5* I < \cn-k\y/g(x) , 
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where 




AlBl 


cos 2 (nT) sin 2 (n7) - 

V , x:=sin 2 (n 7 ) . 


Al 


Bl 


From (10), we obtain |c*| < Ak, k = 1, 2, . . and hence it remains to show that 




A n ~ a n 


By (16), 


bl 


l-< = 1-^* =: /(*) • 


Using standard calculus, one verifies 


< Vg(i) 
/(*) " /(l) 


B n 

— ■ - , 0 < x < 1 

A2 - ’ - - 


and thus (47) holds true, 
b) First, we recall that 


n — 1 


E** = 


fc=i 


X - 1 


, X / 1 


Moreover, for fc = 1, 2, . . ., one has Ak < R k , ak + 6* = r k , and 


„4* 


< 


r 4fc - 1 " r 4 - 1 

Together with (48) (for x = i?/r), we obtain 

n ^ 1 1 n ^ D Ak 

E(f + r )Ak < 4 E<-)‘Tf-r 

— a* ° k k^i r r ~ 1 

Similarly, (49) and (50) (with x = 1/ (J?r)) lead to 


4r 5 


(r 4 — l)(i? — r) v r 


(?) 


2 3 i i 2 3 i 

Br + r)^ < 


.4* 


A:=l 


ajt 6* 


< 


4r 4 


( r ^)* r 4 * “ 1 (r 4 — l)(7?r — 1) 


We prove (46) by verifying that 


n — 1 


a "-* 4 ^ 2r(r 2 + 1) 


k= 1 


(r 2 — !)(/? — r) r 


( 47 ) 


(48) 


(49) 


R n . 
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n ~l 1 

E O n -k . 

r Ak < 

, , a n -kbk 


(~r ■ 


an-kbk ( r 2 - l)(i? - r) V r ; 

The first of these inequalities follows from 

A R k n — Q "~* _ r 2 < ra ~*> + 1 r 2 + 1 , 

* ’ fc> 2 ’ r 2 (n-t) _ 1 - r 2_ x * *-l,...,n-l , 

and ( 48 ) (with x = R/ r). The second one is obtained by making use of 

A t <R k , hz* <x , _L < 2r * , i = , 

On-t 6jt (r 2 - l)r* ’ ’ 

and again ( 48 ) (with x = R/ r). This concludes the proof of the lemma. ■ 

Next, we turn to the derivation of a lower bound for the numbers <r*, / = 1 , . . . , 2 n. 

Using the fact that, by ( 16 ), 

| cos(n7)| / 1 - sin 2 (717) 


< 1 , 


|cos(w7)| _ / 1 - sin 2 ( n7 ) 

yfi \l-sin 2 (n'y )6 2 /B 2 

and part a) of Lemma 2, one obtains from ( 45 ) and (43) the following inequalities: 


*r>l- -^-«n£(|A*(l)| + Ml)|) 


1 an 2 ^n / a n / B n 1 1 w 

6 2 „+^ a n V 

Bnd n , a n £ n |d n |^i a n _ fc , 6„_jt x 


+ i±m + ^- t+ 


* — ** I ** 1 \ / —11 — ^ ''Tl — K v . 

A n b n J“ akb n -k a n -kbk * 


/ = 1 , . . . , 2n 


We set 7 (r/ R) n . With (8) and ( 11 ), one easily verifies that 

b n ^ a n 

B < 7 ’ A < ^ 7 ’ 

/in 

and, together with ( 18 ), the estimates 

’ - 4 n i> n < fl" ’ A„t„ < 27 (51) 

follow. Furthermore, from now on it is assumed that 7 < 1/2, and then 

< r-^TT ( 52 ) 

<*n >ln 1 - 47 2 V 

is guaranteed. By using ( 51 ), ( 52 ), and the inequalities stated in Lemma 2b), we finally 
deduce from ( 50 ) the following lower bound: 

<r,* > 1 — 27 — — (1 + ( 2r2 + 1 )( r2 + j) + 2 ( R ~ r ) 2^\ 

' 7 I-47 2 (r*-l)(R-r)\ r* 7 + r(Rr - l) 7 ) ’ 

/ = 1, 2, . . . , 2n 

In view of Theorem 3 and Theorem 1, this estimate leads to part a) of the following 
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Theorem 4. Let n > 2, c 6 0£r, and R > r > 1 . Then: 

a) q n (z‘, c) is the optimal polynomial for (1 ) with corresponding minimal norm 


D n (r,c) 


r n + l/r" 

.R n + l/.R n ’ 


if 

7 = 7 (r, R, n) := (-^) n 

is such that 7 < 1 / 2 and 


(l- 2 7 ) 2 (l+ 2 7 )> 


8r 5 

(r 4 - 1 )(R - r ) 


( 1 + (JdiM rl+i ) 7 + 


2jR-r) 
r(Rr — 1 ) 



( 53 ) 


holds. 

b) There exists a number Ro(n, r) such that g n (z; c) is the extremal polynomial of (1 ) for 
all 


cEOSr with R>Ro(n,r ) 


c) Let c £ d£n be such that 

9 r 4 - 1 

« > r . ( 54 ) 

Then, there exists an integer «o(r, R) such that q„(z;c) is the extremal polynomial of (1) 
for all n > no(r, R). 

Proof. Only parts b) and c) remain to be proved. For fixed r and n, 7 (r,R,n) — »• 0 if 
R — » 00, and ( 53 ) is clearly satisfied if R is sufficiently large. Similarly, if r and R are 
fixed, the condition ( 54 ) guarantees that ( 53 ) is true if n is large enough. This concludes 
the proof of Theorem 4 . ■ 

Remarks 6. Supported by numerical tests, we conjecture that part c) of Theorem 4 is 
true for arbitrary R > r > 1. 

7 . It follows from R > r > 1 that 


(2r 2 + l)(r 2 + 1) 2(R-r) 

r 4 1 r{Rr-\y 


< 1 + 67 + 27 2 < 


4(1+27) 


for all0<7<l/2. Thus ( 53 ) is true if 7 < 1/ 2 satisfies the stronger condition 


(1 - 2 7 ) 2 > 


18 r 5 

(r 4 - I)(f2 — r) 


( 55 ) 
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Using (55), one easily obtains explicit formulas for numbers R 0 (n,r) with the property 
stated in Theorem 4b). E.g. set 

Ro(n,r) := r max -1 j (56) 

Then, for all R > R 0 (n, r) 


0-2(^)”) J > j > 


18r 5 

(r 4 - 1)(R - r) 


and, in particular (r/ R) n < 1/2. Hence, R 0 (n,r) is suitable for Theorem 4b). 

8. Let G n (r) denote the set of all points c G <F\ £ r for which q„(z;c) is the optimal 
polynomial for (1). By Theorem 4b, G n { r ) is an unbounded set. More precisely, we proved 
that 


C G Gn(r) for all c G (T with |c| > ~(R 0 + — ) , 

2 R 0 

where R 0 = R 0 (n,r ) is given by (56). The boundary of (7 n (r) is a closed Jordan curve 
which, in view of Theorem 3, is composed of pieces of d£ r and of pieces of the curves 


a t( c ) — 0 , l = 1 , 2 , . . . , 2n . 


We have computed these curves numerically for a number of cases. Some typical pictures 
(for r = 2, n = 2, . . . , 5) are shown in Fig. a-d. Because of the symmetry with respect 
to the real and imaginary axis, we have only plotted the first quadrant. G n ( r ) is the region 
exterior to these curves including the parts of its boundary which are described by the 
curves crf(c) = 0. Note that near the real axis the boundary of G n {r) is given by d£ r . 
This led us to the following conjecture: For all n 6 IN and r > 1 there exists a number 
Po ( n , r ) > 0 such that q n (z;c) is the optimal polynomial of (1) for all c G <T \ £ r with 
|argc| < po ( n , r ). 
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